|hep-ph/9707395 
FTUV/97-33 
IFIC/97-32 
July 1997 



00 



X 



Charged Higgs Boson and Stau Phenomenology in 
the Simplest R— Parity Breaking Model 



X5 ■ A. Akeroyd, Marco A. Diaz, J. Ferrandis, M. A. Garcia- Jareno 

[T I I and 

^ ■ Jose W. F. Valle 

(N ■ 



Departamento de Fisica Teorica, IFIC-CSIC, Universidad de Valencia 



(N 
> 

■ Burjassot, Valencia 46100, Spain 

cn 

o 

^ ■ Abstract 

' , We consider the charged scalar boson phenomenology in the simplest effective 

D . low-energy R-parity breaking model characterized by a bilinear violation of R-parity 



in the superpotential. This induces a mixing between staus and the charged Higgs 
boson. We show that the charged Higgs boson mass can be lower than expected in 
the MSSM, even before including radiative corrections. We also study the charged 
scalar boson decay branching ratios and show that the R-parity violating decay rates 
can be comparable or even bigger than the R-parity conserving ones. Moreover, if 
the stau is the LSP it will have only decays into standard model fermions. These 
features could have important implications for charged supersymmetric scalar boson 
searches at future accelerators. 



1 Introduction 



A lot of emphasis has been put into the phenomenological study of the supersymmetric 
Higgs boson sector [|1|]. However, so far most of these phenomenological studies have 
been made in the framework of the Minimal Supersymmetric Standard Model (MSSM) 
0, H with conserved R-parity R-parity is a discrete symmetry assigned as Rp = 
(^—l'^{3B+L+2S) ^ where L is the lepton number, B is the baryon number and S is the spin 
of the state. If R-parity is conserved all supersymmetric particles must always be pair- 
produced, while the lightest of them must be stable. In particular, supersymmetric Higgs 
scalar bosons must decay only to normal standard model particles or to pairs of lowest- 
lying supersymmetric particles, which are usually heavy. On the other hand staus decay 
only to supersymmetric states, like a neutralino and a tau lepton. 

The study of alternative supersymmetric scenarios where the effective low energy 
theory violates R-parity ||5| has recently received a lot of attention p both due to its 
phenomenological interest, as well as due to the intrinsic importance of investigating the 
issue of R-parity breaking at a deeper level. 

It is well-known that the simplest supersymmetric extension of the Standard Model 
violates R-parity through a set of cubic superpotential terms involving a very large number 
of arbitrary Yukawa couplings. Although highly constrained by proton stability, many of 
such scenarios could still be viable. Nevertheless their systematic study at a phenomeno- 
logical level is hardly possible, due to the enormous number of parameters present, in 
addition to those of the MSSM. 

As with other fundamental symmetries, it could well be that R-parity is a symmetry 
at the Lagrangian level but is broken by the ground state. In order to comply with LEP 
precision measurements of the invisible Z decay width these models require the introduc- 
tion of SU{2) (g) U{1) singlet superfields 0. Such scenarios provide a very systematic 
way to include R parity violating effects, automatically consistent with low energy baryon 
number conservation. They have many added virtues, such as the possibility of providing 
a dynamical origin for the breaking of R-parity, through radiative corrections, similar to 
the electroweak symmetry P]. The simplest truncated version of such a model, in which 
the violation of R-parity is effectively parametrized by a bilinear superpotential term 
eiL'^Hl has been widely discussed |^]. It has also been shown recently that this 
model is consistent with minimal N=l supergravity unification with radiative breaking of 
the electroweak symmetry and universal scalar and gaugino masses. This one-parameter 
extension of the MSSM-SUGRA model therefore provides the simplest reference model 
for the breaking of R-parity and constitutes a consistent truncation of the complete dy- 
namical models with spontaneous R-parity breaking proposed previously |]^. In this case 
there is no physical Goldstone boson, the Majoron, associated to the spontaneous break- 
ing of R-parity, since in this effective truncated model the superfield content is exactly 
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the standard one of the MSSM. Formulated as an effective theory at the weak scale, the 
model contains only two new parameters in addition to those of the MSSM. Therefore our 
model provides also the simplest parametrization of R-parity breaking effects. In contrast 
to models with tri-linear R-parity breaking couplings, it leads to a very restrictive and 
systematic pattern of R-parity violating interactions, which can be taken as a reference 
model. 

In this paper we focus on the phenomenology of the charged scalar boson sector 
of the simplest R-parity breaking model. This complements a previous study of the 



electrically neutral sector |TT[. We show that 



1. the mass of the charged Higgs boson can be lower than expected in the MSSM, even 
before including radiative corrections, 

2. if the stau is the LSP it will have only R-parity violating decay channels into stan- 
dard model fermions, 

3. the branching ratio for the R-parity violating charged Higgs boson decays can be 
comparable or even bigger than the R-parity conserving ones. 



We illustrate how these features arising from the charged scalar boson sector of the sim- 
plest R-parity breaking model could play an important role in designing charged super- 
symmetric scalar boson searches at e^e~ colliders such as LEP II. For example they can 
give rise to striking signatures consisting of high multiplicity events, such as di-tau -|- 4 
jets -|- missing energy or 4 taus + 4 jets. Such processes, forbidden in the MSSM, are 
expected to have high rates and negligible background. As for hadron colliders, we also 
can have very high leptonic multiplicity events such as six leptons of which at least two 
are taus, plus missing momentum. This should be easy to see at the LHC, due again to 
the negligible standard model background. 



2 The Model 

The supersymmetric Lagrangian is specified by the superpotential W given by [] 



where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU{2) indices, and e is a com- 
pletely antisymmetric 2x2 matrix, with 812 = 1. The symbol "hat" over each letter 
indicates a superfield, with Qi, Li, Hi, and H2 being SU{2) doublets with hypercharges 
|, —1, —1, and 1 respectively, and U , D, and R being SU{2) singlets with hypercharges 
— |, |, and 2 respectively. The couplings hu, hr, and He are 3x3 Yukawa matrices, and 



*We are using here the notation of refs. H and 111 
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fi and ej are parameters with units of mass. The last term in eq. (HD is the only i?-parity 
violating term. 

Supersymmetry breaking is parametrized with a set of soft supersymmetry breaking 
terms which do not introduce quadratic divergences to the unrenormalized theory |ir 



Vsoft = MfQTQ'; + MfU*Uj + MfD*Dj + Mfl^L'; + MfR*Rj 

+m%^H^*H^ + ml^H^*H^ - [|M,A,A, + ^MXX + ^M'X'X' + h.c] (2) 
+eab [j^ih'^Q'lUjHl + A%h%Q'^^H'^ + A%h%L\R,H'^ - B^iH'^Hl + B^e.^Hl 

and again, the last term in eq. (0) is the only R-parity violating term. 

Following previous discussions we will focus for simplicity on the case of one gener- 
ation, namely the third |p!T|, WM. In contrast we will keep in our discussion the theory as 



defined at low energies by the most general set of soft-breaking masses, tri-linear and bilin- 
ear soft-breaking parameters, gaugino masses and the Higgs superfield mixing parameter 
/i. 

The electroweak symmetry is broken when the two Higgs doublets Hi and if2, and 
the third component of the left slepton doublet Z3 acquire vacuum expectation values. 
We introduce the notation: 

Note that the W boson acquires a mass = \g'^v'^^ where = vl + vl + v1 = 
(246 GeV)^. We introduce the following notation in spherical coordinates for the vacuum 
expectation values (VEV): 

Vi = V sin 9 cos /? 
V2 = V sin 6 sin (3 

V3 = vcos6 (4) 

which preserves the MSSM definition tan/5 = V2/V1. In the MSSM limit, where 63 = w 3 = 
0, the angle 9 is equal to 7r/2. 

In addition to the above MSSM parameters, our model contains three new param- 
eters, €3, V3 and B2, of which only two are independent, and these may be chosen as 63 
and Vs. 



The full scalar potential at tree level is 



Vtotal — ^ 



dW ^ 



+ Vd + Vsoft (5) 
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where Zi is any one of the scalar fields in the superpotential, Vd are the D-terms, and 
Vsoft the SUSY soft breaking terms given in eq. (H). 

The scalar potential contains linear terms 

Vunear=t\x\ + tlxl+tli'^. (6) 

where 

^1 = + ^J^)vi - B - ^le^v^ + \{g'^ + g''^)vi{vl - vl + vl) , 

tl = {ml,+fx^ + el)v2-Bfxvi + B2e3V3-l{g^ + g'^)v2{vl-vl + vl), (7) 

^3 = i^l^ + 4>3 - fJ^e^vi + B2e3V2 + lig^ + g'^)v3{vl - vi + vj) . 

These i = 1,2,3 are the tree level tadpoles, and are equal to zero at the minimum of 
the potential. 

Now a few theoretical comments on the model. The first refers to the choice of basis 
in the original Lagrangian in eq. (|l|). We could have used a rotated basis in which the 
bi-linear coupling disappears [^]. However, if we choose to do so, new R-parity violating 



terms appear not only in the tri-linear superpotential, in the form of a DQL term, but 
also in the scalar sector of the theory, i.e. the Higgs potential due to supersymmetry 
breaking. Authors who adopt this basis have a tendency to neglect SUSY breaking 
in the Higgs potential, which is most crucial for our subsequent analysis. Therefore we 
prefer to use in our calculations the basis in which the bi-linear term is not rotated away. 
The final physics results are completely basis-independent IT] 



Another important feature of this model is that lepton number is violated by the 
es term and by the presence of the sneutrino vacuum expectation value ^3. This induces 
a mass for the tau neutrino since Ur mixes with the neutralinos (see the Appendix). 
This mass turns out to be proportional to an effective neutralino mixing parameter ^ = 
(/if 3+63^1)^ characterizing the violation of R-parity, either through f 3 or 63. One can show 
that this parameter corresponds to the R-parity violating VEV in the rotated basis. If 
we stick to the simplest unified supergravity version of the model with universal boundary 
conditions for the soft breaking parameters, then ^ will be small since contributions arising 
from gaugino mixing will cancel, to a large extent, those from Higgsino mixing. This 
cancellation will happen automatically so that in this case m,^^ will be naturally small 
and radiatively calculable in terms of the bottom Yukawa coupling hh, thus accounting 
naturally for the smallness of the Ut mass in this model. In the language of 



one can 



say that universality of the soft breaking terms implies an approximate and radiatively 
calculable alignment and as a result a suppression in ^ and on m^^. This offers a hybrid 
scenario combining the see-saw and radiative schemes of u mass generation. The role of 
the right-handed mass which appears in the see-saw model is played by the neutralino 
mass (which lies at the SUSY scale) while the role of the seesaw-scheme Dirac mass is 
played by the effective neutralino mixing ^ which is induced radiatively. The Ur mass 
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Figure 1: Tau neutrino mass versus de R-parity violating parameter ^ 

induced this way is directly correlated with the magnitude of the effective parameter ^ 
so that R-parity violation acts as the origin for neutrino mass. In Fig. |] we display the 
allowed values of m^,^. One can see that mj,^ values can cover a very wide range, from 



MeV values comparable to the present LEP limit |]T9[ down to values in the eV range, 
even though the individual values of and £3 can be rather large, e. g. f 3, €3 ~ 100 
GeV, because they are naturally "aligned" , without any need for fine-tuning to get small 
m^^. The alignment follows from the universality of the soft breaking parameters at the 
unification scale [jlO[. However, for generality, in this paper we take the parameters at the 
weak scale independent {i.e. no universality assumption), and always impose m^^ ^ 20 
MeV which corresponds to the laboratory limit. 

This brings us to a discussion about cosmology. Clearly our model leads to a tau- 
neutrino which can be much heavier than the limits that follow from the cosmological 



critical density as well as primordial nucleosynthesis would allow However, in this 
model the u^- is unstable and decays via neutral current into three lighter neutrinos pl |. 
In order for this mode to be efficient we estimate that rriu^ must exceed 100 keV or so. On 
the other hand, to avoid problems with primordial nucleosynthesis it is safer to consider 
masses below 1 MeV or so. In order to sharpen these estimates (which are not strict 
bounds) a detailed investigation is required [^]. 

One should bear in mind, however, that rriu^ can be as large as the present labora- 



tory bound [|T9] in the more complete versions of the model in which R-parity is broken 
spontaneously due to sneutrino expectation values |0, H. This is so because such models 
contain a majoron, denoted as J, which opens new decay channels Uj. into u + J where 
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V 



is a lighter neutrino as well as new annihilation channels + z/,- ^ J + J. It has 



been shown explicitly that the lifetimes that can be achieved in the spontaneous broken 
R-parity versions of the model can be sufficiently short to obey the critical density limit 
Moreover, it has been shown that the annihilation channel is efficient enough in 



order to comply the primordial nucleosynthesis bound |2^, while decays may also play 
an important role 



Finally, a word about the magnitude of R-parity violation. It will depend to some 
extent on the process considered. Some R-parity violating observables turn out to be 
proportional to an effective neutralino mixing parameter ^ = (/if 3 + e^Vi Y characterizing 
the violation of R-parity, either through V3 or €3. An example is the mass of the tau- 
neutrino (see below). However, not all R-parity- violating processes are determined by 
some single production processes or R-parity violating decay branching ratios, such as the 
ones discussed in the present paper, depend separately on V3 or €3 and can be rather large 
even for small ^ and rrii,^ . An obvious and important example is the decay of the lightest 
neutralino, which is determined by 63 only since, in the rotated basis it is determined by 
the D Q L superpotential term only. A detailed discussion lies outside the scope of this 
paper O . 



3 Scalar Mass Matrices 



The mass matrix of the charged scalar sector follows from the quadratic terms in the 
scalar potential 



^quadratic — \Hl i H2 ^^L ^^r] 



'R 



+ ... 



For convenience reasons we will divide this 4x4 matrix into 2x2 blocks in the following 
way: 

(9) 



Ml, M% 



where the charged Higgs block is 



M 



HH 



(10) 



Vl 



and hj. is the tau Yukawa coupling. This matrix reduces to the usual charged Higgs mass 
given by 



matrix in the MSSM when we set vs = = and we call = Bfi. The stau block is 



M 



2 

TT 



6 



1^2 



f3 



'1)3 



*3 

t)3 



We recover the usual stau mass matrix again by replacing = 63 = 0, nevertheless, we 
need to replace the expression of the third tadpole in eq. (0) before taking the limit. The 
mixing between the charged Higgs sector and the stau sector is given by the following 
2x2 block: 

/iea - ^h^vivs + ^g^viv^ -^263 + lg'^V2V3 

and as expected, this mixing vanishes in the limit f 3 = 63 = 0. The charged scalar mass 
matrix in eq. (|^), after setting ti = t2 = ts = 0, has determinant equal to zero since one 
of the eigenvectors corresponds to the charged Goldstone boson with zero eigenvalue. 

For completeness, we give the neutral Higgs sector mass matrices. The quadratic 
scalar potential includes 



Ml, 



(12) 



^quadratic 2 



+ 



(13) 



where the CP-odd neutral scalar mass matrix is 



M 



po 



-Ate3 



V3 _|_ 
Vl Vl 



t2 



Ba^ - ^263^ + 

1^ V2 V2 V2 

—B2t^ 



Vl 



-Ate3 
—-8263 
- -8263 



V2 



V-i 



(14) 



'113 ^ -^iis 

This matrix also has a vanishing determinant after the tadpoles are set to zero, and the 
zero eigenvalue corresponds to the mass of the neutral Goldstone boson. The usual MSSM 
mass matrix of the pseudoscalar Higgs sector is recovered in the limit ^3 = 63 = in the 
upper-left 2x2 block, and the third component corresponding to the imaginary part of 
the sneutrino decouples from it. Note that in this limit, the MSSM mass of the sneutrino 
is recovered provided we replace the expression for the tadpole in eq. (0) before taking 
the limit. 



The neutral CP-even scalar sector mass matrix in eq. (p!3| ) is given by | 

2 _ 
5.0 - 



(15) 



-Bf^ - \glv,V2 + i^2^2 _ 5^^^^ + I B^^^ 

-f^es + \glv,v, 5263 - \glv2v, /xe3^ - B2e, f^ + {gWz + | 

where we have defined (^f = g'^ + g''^. In the upper-left 2x2 block, in the limit f 3 = 63 = 0, 

the reader can recognize the MSSM mass matrix corresponding to the CP-even neutral 

^Note that in the non-diagonal entries of this matrix the terms involving gauge couplings correct from 
those given in ref. plj by a factor 2. 



\9lv2V3 
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Higgs sector. Similar to the previous case, in this hmit the third component decouples from 
the other two and corresponds to the real part of the sneutrino, which become degenerate 
with the imaginary part of the sneutrino. Another way of looking at the separation of 
the sneutrino field into real and imaginary parts is through a sneutrino-anti-sneutrino 45 



degrees mixing ^7 



In the general case, there will be a mixing between the Higgs sector and the stau 
sector. The three mass matrices in eqs. (|^), ([I^) , and (|15D are diagonalized by rotation 
matrices which define the eigenvectors 



ft 



R 



Ht 
Hi 

ft 



T 



R 



"GO" 




>?" 




' h ' 






A 


= Rpo 




1 


H 


= Rs<i 

























(16) 



and the eigenvalues diag(0, my±, m?±, m?±) = .Rs±Mg±i2s± for the charged scalar sector, 

"""i '''2 

diag(0, m^, m?i ) = Ry><'My,„Ry,(, for the CP-odd neutral scalar sector, and diag(m^, m^, m?R) 
RsoMsoRso for the CP -even neutral scalar sector. 

The labelling for the different eigenstates is as follows. In the charged scalar sector 
the eigenstate with zero mass is denoted G^. Among the other three, we call staus the 
two eigenstates Sf with the biggest stau component calculated with {Rf±)'^ + (-R^i)^, 
and by convention m -± < m -± . The remaining eigenstate is called charged Higgs . In 
the neutral CP-odd sector, is the eigenstate P° with zero mass. Among the other two 
eigenstates, the one with largest stau component (i?po)^ is called i>^, and the remaining 
state is the CP-odd Higgs A. Similarly, in the neutral CP-even sector, the state 5° with 
the largest stau component (-R^o)^ is called sneutrino z/^. The other two are the neutral 
Higgs bosons h and H, with ruh < rriH- With this notation is the field which is 
mostly the MSSM charged Higgs, but with a small component of stau, and similarly for 
the neutral Higgs bosons. 

If a 3 X 3 matrix M has a zero eigenvalue, then the other two eigenvalues satisfy 

(17) 



1 1 
-TrM±- 

2 2 



(TrM)^ 



4(MnM22 - + MnM33 - Mfg + M22M33 - MI3) 



The CP-odd neutral scalar mass matrix eq. (0) has a zero determinant, so that its 
eigenvalues m?^ and m\ can be calculated exactly with the previous formula. The same 
can be done with the charged scalar mass matrix in the limit hr = 0. In this case, 
the right stau decouples, and the eigenvalues m'jr± and m?± can be also calculated with 
eq. {\n\). Note that this limit is taken here only for the sake of illustration. In all numerical 
calculations we have used the realistic value for which is fixed through an exact tree- 
level relation eq. (^) given in the Appendix. 

One can determine in the tree-level approximation the minimum of the scalar po- 
tential by imposing the condition of vanishing tadpoles in eq. (^. One-loop corrections 
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change these equations to 

U = t^-SU + T,iQ) (18) 

where ti, with i = 1,2,3, are the renormahzed tadpoles, are the tree level tadpoles 
given in eq. (^, 6ti are the tadpole counter-terms, and Tj(Q) are the sum of all one-loop 
contributions to the corresponding one-point functions with zero external momentum. 
The contribution from quarks and squarks to these tadpoles in our model can be found in 



ref. [|10|- In an on shell scheme we identify the tree level tadpoles with the renormahzed 
ones. Therefore, to find the correct minima we use eq. (|^ unchanged, where now all the 
parameters are understood to be renormahzed quantities. 

We have used the diagrammatic (tadpole) method. Although equivalent to the 



effective potential method for minimization purposes the diagrammatic method is 
better than the effective potential when it comes to calculating the one-loop corrected 
scalar masses [^. Following ref. |^ (see also ref. [^), we work in an on-shell scheme 



where by definition the tree level CP-odd Higgs mass and the tree level IV-boson mass 
mw correspond to the respective pole masses. In this case, the renormahzed charged Higgs 
mass mH± is equal to the tree level charged Higgs mass calculated in the previous section, 
plus the following one-loop contributions: 



Am^± = Re 



Ah+h- im\ + ^w) - AAAifn\) - Awwi^w) (19) 



where Ass{v^), with S = H^, A, W, are self energies. Each self energy is infinite by itself, 
but Am^± is finite. 

Before turning to the numerical study of the charged Higgs mass spectrum, let us 
mention that, throughout this paper, except the case where the parameters have been 
fixed (Fig. 5 and 6), we have taken the MSSM parameters varying in the range : 



0.5 < 


tan (3 


< 90 


OGeV < 


M,M' 


< 1000 GeV 


OGeV < 




< 300 GeV 


-500 GeV < 




< 500 GeV 


-200 GeV < 




< OGeV 



(20) 



and the two R-Parity violating parameters varying as: 

- 200 GeV < es < 200 GeV 
-90 GeV < f3 < 90 GeV (21) 

Note that 'm'jj_^, m'jj^, and B2 are fixed through the tadpole equations given in eq. (|^. 
No big differences are observed if we take /i > , and the sign of B is equal to the sign 
of fi because at the weak scale we have m\ oc ^B. We are interested in relatively light 
charged Higgs , so we take |yu| , \B\ < 200 GeV . Similarly we are interested in relatively 
light staus , and that is why we take niR^ , rriL^ < 300 GeV. 
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^ ^ (o): one loop 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

25 50 75 100 125 150 

niA (GeV) 

Figure 2: Tree level and one-loop charged Higgs boson mass as a function of the CP-odd 
Higgs mass m^. The variation of parameters in the scan is indicated in the text. The 
horizontal dashed line corresponds to the VT-boson mass. 

We now turn to the numerical study of the lowest-lying charged scalar boson mass. 
Our results are illustrated in Fig. ^. In Fig. ^ we display allowed values of charged 
Higgs boson mass as a function of the CP-odd neutral Higgs boson mass rriA- Here the 
main point to note is that m/^± can be lower than expected in the MSSM, even before 
including radiative corrections. This is due to negative contributions arising from the 
R-parity violating stau-Higgs mixing, controlled by the parameter 63. We have varied the 
relevant model parameters in the range given by eq. (pOf ) and eq. (plf ). 

The one loop correction Am^± in eq. ([19|) depends on the soft squark masses Mq, 
Mjj and Md, and the tri-linear soft masses At and Af, . Since they appear only through 
radiative corrections , for simplicity we have taken them degenerate at the weak scale : 
Mq = Mu = Md = 1 TeV and At = A^ = A^. 

An alternative way to display the influence of €3 parameter on the charged Higgs 
boson mass can be seen in Fig. |^. In this figure the curves corresponding to different 
£3 and f3 values delimit the minimum theoretically allowed charged Higgs boson mass 
corresponding to those specific values. These curves are found in a scan where the MSSM 
parameters are varied according to eq. ( ^0]) and the R-parity violating parameters €3 and 
V3 are varied according to the label in Fig. |[ Below each curve no points are found. 
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Figure 3: Minimum charged Higgs boson mass versus tan (3. Each curve corresponds to a 
different range of variation of the R-parity violating parameters 63 and 

The radiatively corrected MSSM prediction is recovered, as expected, when 63 = = V3. 
This is indicated by the sohd hne in the figure. For larger values of the R-parity violating 
parameters one sees that the mass can be substantially lower than the MSSM expectation. 
This is due, again, to the negative contributions arising from the R-parity violating stau- 
Higgs mixing, as discussed above. 



4 Production and Decays of Charged Scalars 

Charged scalar pair production cross section can be calculated with the aid of the ZS^ Sj 
Feynman rule, which is equal to iX^g+g- {p + p')^ where p and p' are the momenta in the 
direction of the positive electric charge flow. The A couplings are equal to \zs+s- — 
Rs±X'zs+s--^^± where the couplings in the unrotated basis are 
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Figure 4: Total production cross section of a pair of (a) charged Higgs bosons and (b) 
light staus, as a function of their mass. The centre of mass energy is 192 GeV. 



The differential cross section for the production of two charged scalars is 

da , , „ , „ X 1 



d cos 9r 



[e e 



S,- Sj 



32ns 



A^/^(l, m'i±/s, ■m'i± /s) sin^ 6c 



(23) 



X 



A z,C\Y S 



+ 



9^ 



rrir 



e2\ \ij2 

zs+s- 



m 



2 \2 



where gy 



2s^r and 9% = \ and A (a, 6, c) = + 6^ + — 2ah — 2ac — 26c. 



In Fig. ^ we plot the production cross section of a pair of charged scalars as a 
function of its mass, considering s/s = 192 GeV. In Fig. ^ we have o"(e~*"e~ — > H^H~) 
and most of the points fall on the MSSM curve. The points that deviate from the main 
curve are due to mixing between charged Higgs and right stau. In fact, if the right stau 
were decoupled from the rest of the charged scalars, the charged Higgs pair production 
cross section would be identical to the MSSM for any value of es or v^, and the reason is 
that the upper-left 3x3 relevant sub-matrix of in eq. (|2^) is proportional to the 

identity. 

In Fig. ^ we have cr(e^e~ — > ^i^^f) as a function of m-±. The points concentrate 
around the two MSSM curves corresponding to the cases where is mainly left stau 
(upper curve) and where is mainly right stau (lower curve). The smallness of the 
right-handed stau cross section relative to the left-handed one is understandable because 
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Figure 5: Mixed H rf production cross section as a function of mH± + rn-T at 192 GeV 
centre-of-mass energy. 

the right-handed stau is an SU{2)l singlet. Again, the points which deviate from these 
curves are due to the mixing between the left and right staus. 

An interesting characteristic of our model is the mixed production e^e" — > H^fi 
which is absent in the MSSM. In Fig. |^ we plot the total mixed production cross section, 
defined by a{e^e~ — H^f^) = a{e^e^ — > H^t^) + o"(e+e^ — > H^fi), as a function 
of the sum of the final product masses + m-± for = 192 GeV. This mixed 

production cross section can be sizable, with a maximum value of the order of 0.12 pb. 

As we have already seen, our model allows strong charged Higgs stau mixing, and 
this can substantially affect both the masses and the couplings. As a result the decay 
branching patterns of the charged scalar bosons can be significantly affected. 

We now turn to a discussion of the charged scalar boson decays. Our first result 
here relates to the stau. In Fig. ^ we display the stau decay branching ratios below and 
past the neutralino threshold. In this figure we have fixed the parameters as 

/i = 370 GeV tan/? = 2 

= -4.7 GeV M = 2M' = 170 GeV 
5 = 40 GeV £3 = 10 GeV 
Ar = 500 GeV tur.^ = 400 GeV (24) 
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Figure 6: Stau branching ratios possible in our model for a particular choice of parameters. 
Note the neutralino threshold below which only R-parity violating decays are present. 

in such a way as to ensure that the lightest neutralino is about 80 GeV in mass and thus 
may be produced as a decay product of a stau produced at LEP II energies. 

In addition, we have chosen €3 = 10 GeV and W3 = —4.7 GeV {9 = 1.59 rad) to 
demonstrate that we don't need large R-parity violating parameters to obtain sizeable 
effects. Otherwise, the choice of parameters is arbitrary . 

Below the neutralino threshold, the stau is the LSP and will have only R-parity 
violating decays, therefore totally un-suppressed, even for the case of small R-parity 
breaking mixing. The main modes of stau decay in this case are into UrT, cs and cb, 
as clearly seen from Fig. |^. Moreover, one sees that for typical values of the R-parity 
breaking parameters, the stau will decay inside the detector. 

Finally, for the case of the R-parity violating charged Higgs boson decays one can 
see from Fig. |^ that the branching ratios into supersymmetric channels can be comparable 
or even bigger than the R-parity conserving ones, even for relatively small values of e and 

Indeed it is explicitly seen from Fig. |^ that, in the region of small tan j3, the R-parity 
violating Higgs boson decay branching ratios can exceed the conventional ones and may 
reach values close to 100 %, since the R-parity-conserving decay H"^ rz/^ is proportional 
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Figure 7: Charged Higgs branching ratios possible in our model for a particular choice of 
parameters. The R-parity violating decay dominates at low tan (3. 

to tan^ f3 and so is usually dominant for larger tan (3. This figure was obtained for a fixed 
choice of parameters, given as 

/i = -100 GeV rriL^ = 140 GeV 

f3 = 2.66GeV M = 2M' = 100 GeV 
mH± = 93 GeV eg = 4 GeV 

A^ = GeV = 400 GeV (25) 

Note that in eq. (EH) we have chosen /i with the opposite sign of eq. (El). This is done to 



see that the trends in Fig. ^ are not a peculiarity of that particular choice of parameters. 
Again we have taken small R-parity violating parameters: £3 = 4 GeV and V3 = 2.66 
GeV. In order to get a ligth charged Higgs with m//± = 93 GeV we take \fi\ and \B\ small. 

We have also checked that a cosmologically safe 1 MeV Uj. would not modify ap- 
preciably our conclusions. For example, we have verified that Fig. 7 remains unchanged 
if we vary m^^ up to 1 MeV. The reason is that, although typically correlated with rriiy^ 
the R-parity violating branching ratio may be large if the Higgs and stau masses are 
relatively close to each other, even for much smaller m^^ values. 

Another way to see that the dominance of R-parity-violating Higgs boson decays is 
not an accident of the above parameter choice is illustrated in Fig. S. The various curves 
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Figure 8: The curves denote the maximum attainable R-parity- violating charged Higgs 
branching ratio versus tan /5. 



denote the maximum attainable values for the R-parity-violating Higgs boson branching 
ratio B{H^ — > 't^xVi- The parameter space is scanned and the curves represent the 
boundaries above which no points are found. The MSSM parameters are varied according 
to eq. (^) and the range of variation of the R-parity violating parameters es and is 
indicated in the figure. In this way and in absolute generality, we demonstrate that even 
for very small R-parity violating parameters the branching ratio B{H^ — > t^xT) "^^^ be 
close to unity, and that in the region of tan (3^1 the decay — > t^v^- dominates. 



5 Discussion 



We have seen in the last section [see Fig. ^ that if the lightest stau is the LSP it will 
decay only through R-parity-violating interactions, to cs or tp^. As a result it leads to 
decay signatures which are identical to those of the charged Higgs boson in the MSSM. 
However, if it is not the LSP the fi is more likely to have standard R-parity-conserving 
decays such as neutralino plus r, leading to signals that can be drastically different from 
those expected in the MSSM and which would arise from x° — > v^Z* or x° — > tW* . 
Unless €3 and f 3 are extremely small, the neutralino will decay inside the detector. For 
the case of stau pair production in e+e~ colliders, such as LEP II, this would imply a 
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plethora of new high fermion-multiphcity events (multi-jets and/or multi-leptons) . For 
example, di-tau + 4 jets + missing energy if both neutralinos decay into jets through 
neutral currents, or 4 taus + 4 jets if both neutralinos decay into jets through charged 
currents. Such processes are expected to have high rates and negligible background. As 
for hadron colliders, we also can have very high leptonic multiplicity events such as six 
leptons of which at least two are taus, plus missing momentum. This should be easy to 
see at the LHC, due again to the negligible standard model background. For a recent 
discussion see, for example, ref. 



If the charged Higgs is the lightest charged scalar boson one can distinguish two 
scenarios. If the R-parity-violating parameter 63 is small, the charged Higgs boson mainly 
decays to standard model fermions, thus conserving R-parity. However ( as pointed out in 



ref. |jTT| for the case of neutral Higgs) it is possible to obtain very large branching ratios 
for R-parity-violating charged Higgs boson decays, even for moderate €3 values, as long 
as the mass difference between Higgs and staus is not too large. This happens because 
the R-parity-violating decay rates are governed by gauge strength interactions, whereas 
the R-parity-conserving ones are determined by Yukawa couplings. 

In the opposite case of large €3 ~ m^y we always expect large branching ratios for 
R-parity-violating charged Higgs boson decay modes. As a result, one expects again a 
large number of novel signatures arising from neutralino decays. These would be the same 
large multiplicity events that we mentioned above for stau pair production. 

Contrary to the MSSM, where stau and charged Higgs boson signal topologies are 
different, in our model they can be identical, for suitably chosen parameters. If one of these 
signals is observed, it will be difficult to know if it comes from a charged Higgs or a stau. 
To disentangle them, measurements of masses and decay rates of both particles would be 
required. Alternatively, less experimental information would be required to separate the 
origin of the signals if the particle mass spectrum were predicted from theory, as happens 
in supergravity versions of this model [|l2 . 



Another interesting feature of our model is the mixed production e^e^ — > H^f^ . 
If m^o < m-± then one can produce interesting signatures like di-tau + di-jets -|- missing 
energy. This is obtained when if^ — > r^z/^ and f"^ — > ''"^Xi — ^ r^qqi'T- Although 
the cross section for this case is typically smaller, one can see from Fig. ^ that for many 
choices of parameters it may be non- negligible. 



6 Conclusion 

In summary we have considered the most salient aspects of the phenomenology of the 
charged scalar boson sector in the simplest effective low-energy R-parity breaking model 
characterized by a bilinear violation of R-parity in the superpotential. We have shown 
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that the mass of the charged Higgs boson can be lower than expected in the MSSM, even 
before including radiative corrections. We have also studied the charged scalar boson 
decay branching ratios and show that the R-parity violating decay rates can be comparable 
or even bigger than the R-parity conserving ones. Moreover, if the stau is the LSP it will 
have only decays into standard model fermions, therefore totally un- suppressed. In the 
opposite case where it is heavier than the lightest neutralino one expects a plethora of 
exotic high fermion multiplicity events for which the standard model backgrounds should 
be rather small or absent. A detailed analysis of the detectability prospects of the related 
signatures at future accelerators lies outside the scope of the present paper and will be 
taken up elsewhere. 
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Appendix 

For completeness, in this Appendix we collect the Feynman rules relevant for the study 
of charged scalar decays into two fermions, where these two fermions are a chargino-tau 
and a neutralino-neutrino. We also give the exact formula, eq. (PS]), for the tau-lepton 
Yukawa coupling h^. which differs in this model from that of the MSSM. First we set our 
conventions in the fermionic sector. 

As we have already learned, the charged Higgs states Hf and H2 mix with the stau 
states and fj^ and form a set of four charged scalar eigenstates Sf with i = 1,2, 3, 4. 
The charged scalar mass matrix M|± in eq. (P) is diagonalized by a 4 x 4 rotation matrix 
Rs± defined in eq. (0). 

In a similar way, charginos mix with the tau lepton forming a set of three charge 
fermions F^, i = 1,2, 3. In a basis where ^'^'^ = {—iX~^, H\, r^) and ip~'^ = Hf, r£), 

the charged fermion mass terms in the lagrangian are 




(26) 
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where the chargino/tau mass matrix is given by 

M 

Mc ' 



7^9V2 
-^3 







(27) 



and M is the SU{2) gaugino soft mass. We note that chargino sector decouples from 
the tau sector in the hmit £3 = f 3 = 0. As in the MSSM, the chargino mass matrix is 
diagonahzed by two rotation matrices U and V The tau Yukawa couphng hr is chosen 
such that one of the eigenvalues is equal to the tau mass. This is calculated from the 
vacuum expectation values of the model through an exact tree level relation given by 



f + g{e3,V3) 
f-^h{e3,V3) 



(28) 



where 



-e3V3g''{Mv2 + jJ-Vi) 

h{e3,V3) = esvs {fiv,{M^ - M') + Ig'Mvlv^) + \elvl (M^ - 
+ \vl [m^{M^ - M2) - g'M^ivf + it;2) - g'MfiViV^ 



\g^vlvl 



+ \e3vlg'Mv2 - \vig\M'^ 



In our model, the tau neutrino aquires mass, and this is due to a mixing between the 
neutralino sector and the neutrino-tau, forming a set of five neutral fermion i — 
1, ...5. In the basis '0°'^ = (— iA', — iA^, Hi, the neutral fermions mass terms in the 

lagrangian are given by 



1 



{i^^Mi^i,^ + h.c. 



(29) 



where the neutralino/neutrino mass matrix is 



M' 


-\g'vi 
\g'v2 
-\g'vz 





M 

\gvi 
-\gv2 



-\g'vi 


\g'v2 


-\g'v3 


\gvx 


-\gv2 


\g^^ 















£3 





£3 






(30) 



and M' is the f/(l) gaugino soft mass. This neutralino/neutrino mass matrix is diagonal- 
ized by a 5 X 5 rotation matrix N such that 



N*MnN ^ — diag(m^o, m^o, m^o,m^o,m,^^) 



(31) 



where by definition the eigenstate is the neutrino-tau, i.e., with the largest tau com- 
ponent {Ni^y. 
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Now we are ready to work out the SfF^Fj^ Feynman rules. These vertices are 
denoted by 

S^Fj F° — > i\^^^^_pgPj^ + iX^^^^_pQPj^ (32) 

where Pl = — 75) and Pr = 1(1 + 75) are the usual left and right proyection operators. 
For simplicity, we work in a basis where the chaged scalars are unrotated. In this basis 
S"+ = {Hf*, Hi, f^, f^) the vertices are denoted by 



'^i ^ T"^s+F-FO^L + T"^s+F-FO^R 

The relation between A and A' is given by 



(33) 



\-L ijk -ryil \>L Ijk 

■^S+F-FO ~ "^S^'^S+F-FO ' 



\Rijk -pi; \lRljk 

■^S+F-FO — "^S^'^S+F- 



FO 



(34) 



and each of the A' can be read from the following Feynman rules 



Hl*FrF^ 



HlFrF^ 



ftF-F^ 



f + P- F° 



■ 9 
'2 



;i-75)+^YiV,5\/i3(l + 75) 



V2 



9 



:i + 75) 



(35) 



■ 9 
-I — 
2 



;i-75)-^YiV,3F.3(l+75) 



—I 



^y,3iv,i(i + 75) + {n*,u:^ - n;,u:,) (1 - 75) 



The reader can check that from eq. (|35|) we can recover the MSSM Feynman rules taking 
the appropriate hmits. 
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